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1.  Introduction 

The  advent  of  third-generation  synchrotron  radiation  (SR)  sources,  the  development  of  ultra  large  scale  integrated  (ULSI) 
circuits,  and  the  expanding  development  of  space  observatories  have  created  an  increasing  need  for  high  quality,  well 
characterized  optical  components  and  systems  for  use  in  the  soft  x-ray  (SXR)  region.  For  high-brightness  "third  generation" 
SR  sources,  beamline  optics  are  required  to  withstand  extreme  power  loads  while  providing  precisely  tailored  radiation  to 
experimental  stations.  The  ever-increasing  density  of  integrated  circuit  patterns  is  approaching  fundamental  limits  of  what  can 
be  accomplished  with  ultraviolet  excimer  lasers,  and  next-generation  microelectronics  will  require  the  development  of 
sophisticated  and  complex  SXR  optical  systems.  Next-generation  SXR  space  observatories  also  will  require  utilization  of 
complex  aspheric  optical  elements  in  their  telescopes  and  spectrometers  to  achieve  large  aperture  combined  with  the  highest 
possible  spatial  and  spectral  resolutions.  It  should  also  be  noted  that  the  need  for  new  high  precision  SXR  optics  and  optical 
systems  has  pressing  applications  in  many  other  fields,  ranging  from  biology  to  materials  research. 

Due  to  short  operating  wavelengths,  SXR  optics  must  meet  very  stringent  metrological  requirements  not  usual  to 
conventional  optical  elements  used  at  longer  wavelengths.  The  current  lack  of  sufficiently  precise,  yet  practically  convenient 
metrological  technologies  presents  varying  degrees  of  difficulty  in  judging  the  quality  of  SXR  optical  components  and  the 
accuracy  of  an  assembled  optical  system.  At  this  stage  of  the  development  of  SXR  optics,  it  is  desirable  for  manufacturers  as 
well  as  users  to  have  a  practical  design  method  as  well  as  a  simulation  method  that  can  provide  tolerance  estimates  for  the 
fabrication  of  optical  components  needed  to  secure  the  highest  possible  performance  of  a  planned  SXR  optical  system. 

To  contribute  towards  the  achievement  of  these  ends,  we  have  developed  practical  methods  which  incorporated  analytic 
formulas  and  appropriate  definitions  into  a  ray  tracing  method  based  on  Fermat's  principle.  They  were  found  to  be  effective  not 
only  in  designing  a  multi-element  optical  systems  but  also  in  estimating  tolerances  for  fabrication  and  assembly  errors  of  optics 
and  maximizing  the  performance  of  resulted  systems.  The  principles  of  the  methods  and  some  applications  are  given. 

2.  Fermat's  Principle  and  Ray-Tracing  Formulas 

The  theory  of  geometrical  optics  can  be  developed  on  the 
basis  of  a  single  hypothesis,  known  as  Fermat's  principle  [  1  ]. 

For  a  diffraction  grating,  the  light  path  function  F  for  a  ray 
originating  from  a  point  A  and  diffracted  at  a  point  P  ( E„w,f) 
on  the  nth  groove  toward  an  image  point  B  is  defined  by 
F  -  AP  *  PB  *  nmX.  Here,  m  is  the  spectral  order  and  X  is 
the  wavelength  of  the  ray  PB.  Since  the  surface  of  the  grating 
blank  is  expressed  by  $  F  is  also  a  function  of  w  and 

/.  Application  of  Fermat's  principle  to  the  light  path  function 
yields  6 F  -  (dF/d w)6w  ♦  {dF/dl)  hi  -  0.  Fermat's 
principle  as  applied  to  the  diffraction  grating  is  now  stated  as 
( dF/dw )  -  0  and  ( dF/dl )  -  0 ,  and  more  explicitly  as 

dF/dw  -  (L-L%dl/dw)  *  (M-W)  ♦  mX(dn/dw )  -  0, 

(1) 

dF/dl  -  (L-L%d^/dl)  *  (.N-N')  *  mXCdn/dl)  -  0, 

where  (L,M, N)  and  (L'M'.N')  are  the  direction  cosines  of  the 
incident  ray  AP  and  those  of  the  diffracted  ray  PB  of 
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wavelength  X  in  with  order,  respectively,  and  they  are  defined  as 

L  «  (l-x)/AP,  M  =  ( w-y)/AP ,  N  =  (l-z)/AP,  L'  =  ( x'-l)/PB ,  M‘  =  (y'-w)IPB,  N‘  =  (. z'-l)/PB .  (2) 


Solving  these  equations  simultaneously,  we  obtain 

I'  -  £  ♦  T,  M'  -  M  *  m X—  -  T~, 

dw  dw 


N‘  -  N  *  mX  —  -  A 

31  31 


(3) 


A1 1  the  equations  derived  above  hold  for  gratings  of  any 
type,  and  the  direction  cosines  ( L  'M'.N')  of  the  diffracted 
ray  PB  can  be  computed  whenever  the  surface  figure  L,  = 
j[yv,t)  and  the  groove  pattern  n  =  h(w,l)  are  known. 

We  consider  two  image  planes  E  and  E0,  both  of  which 
are  perpendicular  to  the  x-y  plane  and  pass  through  a  point 
BoC^o.O)  on  the  diffracted  principal  ray  of  wavelength  X  in 
with  order,  the  incident  principal  ray  being  A^O.  The  plane 
E0  is  perpendicular  to  the  diffracted  principal  ray  OB„.  The 
plane  E  makes  an  angle  <}>  with  the  plane  E0: 

*'cos((30  ♦  4>)  +  >''sin(p0  ♦  <j>)  -  r0'cos<j>.  (4) 


Fig.  2.  Schematic  diagram  of  two  image  planes  E  and  E0 


Here,  p0  is  the  angle  of  diffraction  of  the  principal  ray  A0OB0, 

<j>  =  0  for  a  monochromator,  <j>  =  x  ~  P0  _  71  ^2  for  a  flat-field 
spectrograph,  and  a  series  of  image  planes  E  with  (j>  =  po 
represents  the  focal  plane  for  a  Rowland-circle  spectrograph.  Thus,  the  intersection  B(x  \y  \z ")  of  the  ray  PB  with  the  image 
plane  E,  i.e.,  the  coordinates  of  the  ray-traced  spot  are  determined: 


x'  -  5  ♦  L' d,  y’  .w.M' d,  z'  -l  *N'd, 

d  -  [r0'cosc))-  5 cos(P0t<j))  -  wsinf  pot(j))]/[L'cos(  P0*^>)-A/'sin(p0«<t))] . 


(5) 


For  convenience  we  introduce  in  the  plane  E  a  new  rectangular  coordinate  system  whose  Y  axis  lies  in  the  x-y  plane  and 
Z  axis  is  parallel  to  the  z  axis.  The  ray-traced  image  point  B(x  \y  \z  is  then  expressed  in  the  new  coordinates  as  B(Y,Z): 

Y  -  O'- i-osinpo)sec(p0.  4>),  Z  -  z' .  (6) 


3.  Hybrid  Design  Method 


It  is  desirable  to  have  a  practical  method  of  designing  a  multiple-element  optical  system  consisting  of  aspheric  mirrors  and 
varied-line-space  and  curved-groove  gratings,  both  ruled  and  holographic.  For  this  we  have  developed  a  hybrid  design  method 
which  takes  the  advantages  of  both  the  ray- tracing  and  analytic  methods  [2], 


The  hybrid  method  optimizes  design  parameters  by  minimizing  an  analytic  merit  function  Q=  E,  Q(X,)  whose  component 
Q(X)  closely  represents  the  variance  of  the  spots  formed  when  an  infinite  number  of  rays  of  a  design  wavelength  A.,  are  traced 
through  the  optical  system  concerned  [3].  A  component  of  the  merit  function  is  defined  as 


_L_f™(™[™(Y-Ffdw<Udz 

WLH  J-wn  Pm  J-m 


g  rWtl  rill  rHtl  ^2 

WLH  J-wn  J-m  J-un 


Aw  Al  iz , 


(7) 
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where  A,.  is  a  wavelength  chosen  in  the  required  scanning  range,  W  and  L  are  the  width  and  height  of  the  illuminated  portion 
of  the  ruled  area,  respectively,  H  is  the  entrance  slit  height,  and 


J-  fwn  fm  fHn  T  dwdldz , 
WLH  J-wn  J-i/2  J-m 


<lh\)  -  ~  ^[w\2fL  *  9/100/300)  *  2L*f0la  ,  2  !*«/«]  .  J-  [2  -  H2fm)  *  i2tf2/on 


960 


W'[l  W2f300  ( L  2fm  *  H2fm)  .  L*fil 0  ♦  H*&\  .  ^  W&  .  -JL  W'LWaZt  *  2  , 


1728 


(8) 


*  -jj  (i : 2  Sow 'H 1  gm  )  ~  {  W2  L 2  (sL *2  golo  gna)+H2[w2(  g2,  ♦  2  gool  g20l  )  ♦  2  L 2  (  g010  gou  .  ^  g02,  )  ] } 
^L*g0l0g030.^\lV2L2g21o(W2gm.2L2goio),L2H2gon(2L2g0io.H2g012)tHHw*g^tL*g^)} 

*  448  ^  ^030  *  864  ^  ^  H  ^ &10 £012  * £021  £201 )  •  .  '  . 

In  brief,  for  a  given  optical  geometry,  we  (1 )  generate  nine  rays  («  =  1 ,2,-,9)  of  wavelength  A,  randomly  from  the  source, 
(2)  determine  the  intersections  of  nth  ray  with  the  surfaces  of  the  grating  and  the  image  plane,  (w„,/„)  and  (Yn,Z„),  by  means  of 
ray  tracing,  (3)  substitute  the  values  of  (w„./„),  (F„,Z„),  and  the  height  zn  of  the  source  point  into  the  analytic  formulas  for  the 
spot  diagrams  nig-.?  oT--' 


Yh  -  VlOO  +  *V4 0  *  lnftati  '  lnZnfm  ’  Z»/oo2  ♦  wnfl 00  *  W»/n/I20  +  *nhZJ\ til  *  W»*»/102  *  ™. 

ZK  ■  *„So01  *  '.*010  *  ”nln8 110  +  3V«*101  *  U210  ♦  ^201  *  ^OJO  *  4^  +  l£im  +  “>  ! 


(9) 


(4)  solve  the  nine  simultaneous  equations  fork's  and  gijk's,  (5)  repeat  steps  (1 )-(4)  for  the  other  design  wavelengths,  and  (6) 
calculate  <2(A,)’s  with  the  values  of  the  Jjks  and  giJk s.  Finally,  the  design  parameters  are  determined  by  minimizing  the  merit 
function^?  =  E,  Q(A,).  If  necessary,  we  may  use  Q  =  2,  e(A,)Q(A,),  e(A,)  being  a  weighting  factor.  h* 

4.  Resolving  Power 


4. 1.  Spectroscopic  Systems 


visw 

;lb  Vtffitl® 


‘Of  j»  .if: 


The  resolving  power  of  spectroscopic  systems,  9t=A/A  A,  may  be  defined,  in  accordance  with  the  Rayleigh  criterion,  by 
assuming  that  two  similar  spectrum  lines  of  A  and  A+A  A  are  resolved  when  their  contours  cross  at  the  point  of  each  whose 
height  is  4/tc  2  of  the  maximum  (see  Fig.  3a).  In  this  case,  A  A  is  the  full  width  of  the  ray-traced  line  profile  at  its  (4rit2)-maximum 
point.  The  resolving  power  thus  defined  is  denoted  by  SR(  1 ).  This  definition  would  be  appropriate  only  when  spectral  images 
have  symmetric  profiles  of  Gaussian  type.  To  take  into  account  the  effect  of  asymmetric  line  profiles,  we  define  the  resolving 
power  9t(2)  by  introducing  an  effective  Gaussian  line  profile  in  place  of  an  actual  asymmetric  line  profile  (refer  to  Fig-  3b  for 
defining  the  effective  Gaussian  profile)  [2], 


A  A+AA 


A 


imnaottom  p>« 
in  be**  .griivfe 

,  ,  ;A-t-AA 


(a)  (b) 

Fig.  3.  Definitions  of  the  resolving  power  based  (a)  on  the  Rayleigh  criterion  and  (b)  on  the  effective  Gaussian  profiles. 
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To  define  the  resolving  power  9t( 2),  we  (1)  calculate  the  standard  deviation  oy  of  ray-traced  spots  in  the  direction  of 
dispersion,  (2)  convert  a  Y  to  the  standard  deviation  oA  of  the  spectral  spread  for  the  rays  of  wavelength  A  by  multiplying  the 
reciprocal  linear  dispersion  at  A.,  (3)  represent  the  line  profile  under  consideration  by  an  effective  Gaussian  profile  whose 
standard  deviation  is  given  by  oA,  and  (4)  assume  that  two  similar  spectral  lines  of  A  and  A+A  A  are  resolved  when  their  effective 
Gaussian  lines  are  separated  by  A  A  =2. 643  o  ,  to  make  the  minimum  resultant  intensity  between  the  lines  8/it2  as  great  as  the 
resultant  intensity  at  the  central  maximum  of  either  of  the  lines.  Thus  the  resolving  power  3t(2)  is  defined  by 

»(2)  -  A/AA  -  A/2.643  or  (10) 

4.  2.  Demagnifying  Projection  Lithography  Systems 

The  square-wave  modulation  transfer  function  (MTF)  has 
frequently  been  used  to  evaluate  the  imaging  characteristics  of 
demagnifying  projection  optics  [4],  On  referring  to  our  previous 
experience  gained  through  projection  lithographic  processes 
including  resist  processes,  we  assume  that  a  projected  line-and- 
space  pattern  at  a  given  spatial  frequency  is  properly  resolved 
when  the  modulation  of  the  intensity  distribution  of  the  mask 
image  formed  on  a  wafer  is  45%  or  more. 

To  estimate  the  performance  of  a  demagnifying  projection 
optics  by  means  of  ray  tracing,  we  ( 1 )  consider  a  grid  pattern  of  a- 
pm  line-and-space  on  the  reflection  mask  (see  Fig.  4),  (2)  assume 
that  the  space  portion  of  this  grid  pattern  emits  light  into  the 
entrance  pupil  of  the  system,  (3)  choose  three  6a  x  6a- pm2  areas 
centered  around  the  two  comer  points  at  the  right  (or  left)  and  the 
central  point  of  the  ring  field  formed  on  the  mask,  and  take  these 
areas  as  distributed  light  sources,  each  containing  9  elemental 
sources  of  a  x  a  pm2,  (4)  generate  rays  randomly  from  each  one 
of  these  three  areas  and  construct  spot  diagrams  in  the  plane  of  the 
wafer  by  tracing  rays  through  the  projection  optics,  (5)  construct 
the  intensity  distributions  of  the  projected  image  in  the  two 
directions,  x  and  y,  parallel  to  the  grid  lines  by  counting  the 
number  of  spots  falling  into  individual  zones  of  (6aA//100)-pm  wide  taken  in  the  x  and>>  directions,  M  being  the  magnification 
of  the  projection  optics,  and  (6)  estimate  the  modulation  of  the  intensity  profiles  of  the  pattern  thus  obtained  in  the  two  directions 
and  apply  the  "s45%"  criterion  to  the  estimated  modulations. 

5.  Spectral  Purity 

The  spectral  purify  of  the  beam  emerging  from  the  exit  slit  is  an  important  factor  in  evaluating  the  performance  of  a  high- 
resolution  monochromator  on  a  synchrotron  radiation  beamline,  whose  source  is  of  continuum.  It  is  also  important  that  the 
resolving  power  estimated  should  not  contradict  with  the  spectral  purity. 

The  spectral  purity  may  be  examined  by  analyzing  the  frequencies  at  which  rays  of  various  wavelength  components  of  the 
continuum  pass  through  the  exit  slit  of  a  given  width.  Referring  to  Fig.  5,  we  (1)  assume  that  the  monochromator  is  tuned  to 
a  wavelength  A0  and  that  a  continuum  band  which  extends  over  a  range  of  A,<  A  <:  A2  is  incident  on  the  grating,  choice  of  the  limits 
A,  and  A2  being  made  on  account  of  the  reciprocal  linear  dispersion,  (2)  generate  rays  of  various  wavelengths  in  the  continuum 
band  randomly  in  the  bending  magnet  source  or  the  central  cone  of  the  undulator  under  consideration  and  trace  the  generated 
rays  through  the  beamline  optics  until  the  number  of  rays  passed  through  the  exit  slit  reaches  1000,  (3)  identify  each  through-ray 
by  its  wavelength  A  and  the  Y  coordinate  of  its  intersection  with  the  plane  of  the  exit  slit,  the  Y  axis  being  taken  as  the  line 
perpendicular  to  the  slit  and  going  through  the  center  of  the  slit  opening,  (4)  construct  the  spectral  purity  plot  diagram  by  plotting 
(A,T)'s  for  the  through-rays,  (5)  estimate  from  the  spectral  purity  plot  diagram  the  frequency  distribution  v(A)  of  the  wavelength 
component  A  by  counting  the  number  of  points  (A,T)'s  falling  onto  individual  parallel  zones  having  a  width  of  (A2-A,)/100  run 
and  a  height  equal  to  the  width  of  the  exit  slit,  and  (6)  calculate  the  spectral  purity  A/Ar,  where  AR  and  A  T  are  the  area  under  the 
frequency  distribution  curve  between  (A0-1.3215oa)  nm  and  (A0+1.3215oa)  nm  and  the  total  area  under  the  whole  curve, 
respectively. 


a  pm 


< 6  a  pm  


Fig.  4.  Source  pattern  on  the  mask. 
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Fig.  5.  Intensity  distribution  of  the  continuum  band  and  explanations  of  through-rays  and  spectral  purity  plot  diagrams. 

6.  Slope  error  and  figure  error 


We  consider  an  ellipsoidal  mirror  whose  surface  figure  is  defined  by 


l  -  A  -  A 


A« 


(11) 


in  a  rectangular  coordinate  system  attached  to  the  mirror.  The  vertex  O  of  the  mirror  is  taken  as  the  coordinate  origin  and  the 
mirror  normal  at  O  as  the  x  axis.  In  Eq.  (1 1),  A,  B,  and  C  are  the  ellipsoidal  semiaxes  with  respect  to  the  x,  y,  and  z  axes, 
respectively,  and  (t,w >,[)  are  the  coordinates  of  a  point  P  on  the  mirror  surface.  The  quantity  A  §  is  the  deviation  from  the  ideal 
ellipsoidal  surface,  and  is  assumed  to  be  constituted  of  the  surface  figure  error  A  \FE  and  thermal  deformation  A 

A£  -  A£j ^  ♦  A (12) 


The  figure  error  arises  from  the  imperfection  in  polishing  the  surface  of  a  mirror  or  a  grating  blank,  and  it  often  is 
represented  by  values  of  the  slope  error.  We  assume  that  the  slope  errors  sy  and  s,  at  P(£,w’,/)  in  they  and  z  directions  are 
distributed  randomly  with  a  probability  density  function 


A‘y,*z) 


1 

- exp 

2it  aj02 


(13) 


where  o ,  and  o2  are  the  standard  deviations  of  sy  and  sz,  respectively.  We  assumed  here  that  the  sy  and  sz  have  no  mutual 
correlation  and  their  respective  averages  over  the  mirror  surface  are  zero.  The  sequence,  {  (sy,  sz), },  of  normal  random  slope 
errors  with  the  probability  density  J(sysz)  are  generated  from 


0>),  =  [d(A£*E)/dH’].  -  o,(-21ogeH(.)1/2cos  2ji  Uj.j, 
(*z>,  -  p(A  5JB)/a/]<  -  o2(-21og<«I.)1/2sm2Jt«I.,1, 


(14) 


with  the  aid  of  the  sequence  {  m,  }  of  uniform  random  numbers  on  the  interval  [0,1]. 

The  corresponding  figure  error  A  E,fe  is  obtained  by  integrating  d(A£re)  =  s^dir  +  sfdl  under  an  assumption  o,  -  o2  =  oSE. 
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PE  -  oSB(-21og#«,)1/2(wcos2  5iH,,1 


/sin  2  3i  u(tl). 


(15) 


1.  SR  Beamline  Optics 

Examples  are  given  for  the  design,  resolving  power,  spectral  purity  of  an  undulator  beamline  optics,  and  thermal 
deformation. 

7.  1.  Validity  of  various  ray-deviation  formulas 

Equations  (9)  express  the  aberrations  (up  to  third  order)  in  spectral  images  formed  in  the  plane  S  (see  Fig.2).  They 
represent  ray-traced  spot  diagrams  also,  and  they  are  called  spot  diagram  (SD)  formulas.  The  accuracy  of  Eqs.  (9)  is  limited 
only  by  the  neglect  of  the  fifth-  and  higher-order  terms  in  performing  power  series  expansions. 

It  is  of  practical  interest  to  examine,  with  reference  to  exact  ray  tracing,  the  validity  of  the  SD  formulas  and  also  that  of 
similar  formulas  based  on  light  path  function  (LPF).  The  LPF-based  formulas,  which  are  widely  in  use  for  the  design  of  grating 
instruments,  are  given  by  [5] 


Y  -  r0' sec  p 0(dF/dw),  Z  -  -  (zr07r)  *  r^dF/dl), 


(16) 


where  the  image  plane  is  assumed  to  be  perpendicular  to  the  principal  diffracted  ray  (E0  in  Fig.  2  and  <j>  =  0  in  Eq.  (4)).  The 
expansion  coefficients  of  the  light  path  function  F  contain  the  coordinates  (r',  p.z')  of  point  B,  which  are  unknown  functions  of 
the  coordinates  (£,h>,/)  of  point  P  (see  Fig.  1 ).  Partial  differentiation  of  F  therefore  is  not  possible  in  principle.  To  circumvent 
this  problem,  it  is  customarily  assumed  that  the  deviation  of  the  spot  B  formed  by  the  ray  APB  from  the  spot  B0  of  the  principal 
ray  AOB0  in  the  image  plane  is  negligible  so  that  the  unknown  quantities  r\  p,  and  z'  in  the  expansion  coefficients  of  the  light 
path  function  can  be  approximated  by  the  corresponding  quantities  r0',  p0,  -zr07r,  which  are  all  associated  with  the  principal 
ray  A„OB0 : 


F  ~ 


[i?]r'.r#',p-p0,r'-*r07r- 


(17) 


With  this  approximation,  Eqs.  (16)  are  expressed  as 
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It  can  be  shown  analytically  that  the  LPF-based  formulas  are  correct  to  third  order  only  when  P0  =  0  and  the  tangential  and 
sagittal  focal  curves  cross  on  the  grating  normal  [1],  The  same  results  are  obtained  [1]  for  another  ray-deviation  formulas  that 
are  based  on  the  wavefront  aberration  (WFA)  theory  [6,7],  In  practice  these  LPF-  and  WFA-based  formulas  can  safely  be 
applied  to  cases  in  which  defocus  (F200)  and  astigmatism  (f020)  are  small  over  the  design  wavelength  range  (i.e.,the 
corresponding  angles  of  diffraction  P0  are  not  too  large). 

The  validity  of  individual  formulas  is  depicted  clearly  in  Fig.  6,  which  compares  the  spot  diagrams  generated  for  a  142° 
constant-deviation  grazing  incidence  monochromator  by  the  SD-  and  LPF-based  formulas  with  those  constructed  by  exact  ray 
tracing.  Specifications  of  the  monochromator  are:  wavelength  range,  15nm-  150  nm;  m  =  +1;  r—  319.9  mm;  r0'  =  317.8  nrn; 
R  =  1000  mm;  l/o  =  550  grooves/mm;  2 a  =  3.877791  x  10'9  mm;  6b  =  -6.844844  x  10‘14  mm;  4c  =  3.654780  x  10'2°  mm;  and 
ruled  area,  60  (W)  x  30  (L)  mm2.  The  monochromator  has  a  rather  large  astigmatism  as  can  be  seen  in  Fig.  6.  The  spot 
diagrams  in  Figs.  6b  (SD-based  formulas)  are  very  similar  to  the  corresponding  ones  in  Fig.  6a  (ray  tracing),  proving  the  validity 
of  the  SD-based  formulas.  On  the  other  hand,  the  spot  diagrams  generated  from  the  LPF-based  formulas  (Fig.  6c)  are  quite 
different  in  shape  and  magnitude  from  those  constructed  by  exact  ray  tracing  and  the  SD-based  formulas.  Spot  diagrams 
constructed  by  the  WFA-based  formulas  also  deviate  considerably  from  Figs.  6a  and  6b  (see  Ref.  1).  These  results  show  the 
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limitation  in  the  applicability  of  the  LPF-  and  WFA-based  formulas  to  a  system  in  which  astigmatism  is  not  corrected  to  a 
sufficient  degree. 


A=20nm  60nm  lOOnm  I40nm 


IMAGE  WIDTH  (mm) 


Fig.  6.  Spot  diagrams  constructed  for  the  1 42  °  constant-deviation  grazing  incidence  monochromator  with  (a)  exact  ray  tracing, 
(b)  SD  formulas,  and  (c)  LPF  formulas. 


7.  2.  Design,  resolving  power,  and  spectral  purity  of  undulator  beamline  optics 
A.  Design 

We  consider  an  objective 
Monk-Gillieson  (M-G)  type 
grazing  incidence  VLS  plane 
grating  monochromator  (PGM)  _  18  ill 

shown  in  Fig.  7.  A  fixed  g  '  * 

spherical  concave  mirror  M  0  1 

accepts  radiation  from  a  source 

point  S  at  an  angle  of  incidence  G 

88°  and  delivers  a  converging 
beam  onto  a  ruled  VLS  plane 

grating  G  at  an  angle  of  Fig.  7.  Schematic  diagram  of  the  M-G  type 

incidence  a.  The  diffracted  rays 
of  wavelength  A  in  /wth  order  are  focused  onto  the  fixed  exit  slit  Ex  after  reflected  by  a  movable  plane  mirror  M1.  Wavelength 
scanning  is  carried  out  by  combining  simultaneous  translation  and  rotation  of  M’  with  simple  rotation  of  G  so  that  the  direction 
of  the  exiting  beam  remains  unchanged. 

.  'ii  f\  rflfyi*?  I  jut  T'TffX  ,*^'4  ; 

This  monochromator,  M-G  type  VLS  PGM,  was  designed  by  means  of  the  hybrid  method  described  in  Section  3.  The  design 
conditions  and  the  results  are  summarized  below. 

Given  conditions:  r  =  SM  =  18000  mm,  0  =  88°,  R  (radius  of  curvature  of  M)  =  155700  mm,  D  =  MG  =  200  mm,  <j t-K=  1  °, 
o  =  1/2400  mm,  ruled  area  =  60  (W)  x  25  (L)  mm2,  spectral  order  m  =  +1,  and  wavelength  range  =  0.5-10  nm. 
Designed  parameters:  2a  =  1. 762906*  10"10  mm,  66  =  5.299280*  10'17  mm,  4c  =  1. 354672*  10‘23  mm,  D'  =  GM' =  401.10  mm 
at  A=5  nm,  and  r‘  =  M'Ex  =  1647.78  mm  at  A=5  nm.  (Note  here  that  D  and  r'  are  wavelength  dependent.) 

We  also  designed  holographic  plane  gratings  for  the  monochromator  under  the  same  given  conditions  as  described  above 
for  the  ruled  VLS  gratings.  We  designate  the  holographic  gratings  (HG)  recorded  with  (1)  two  spherical  wave  fronts  and  (2) 
an  aspheric  and  a  spherical  wave  front  as  HG(S+S)  and  HG(A+S),  respectively,  and  the  ruled  VLS  gratings  as  RG(VLS).  The 
design  parameters  thus  obtained  with  A0=44 1 .6  nm  are  summarized  below. 
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HG(S+S):  rc  =  -2194.42  mm,  y  =  -71.2090°;  rB  =  959.60  mm,  6  =6.4964°. 

HG(A+S):  i?,  =  1200  nim ,pc  =  585.44  mm,  qc  =  614.96  mm,  tic  =  -35.3665°,  y  =  -36.8560°;  =  979.50  mm,  6 

=  27.3893°. 

The  performance  of  all  the  gratings  thus  designed  was  evaluated  in  terms  of  resolving  power  and  spectral  purity  determined 
from  ray-traced  spot  diagrams,  assuming  Lu  (undulator  length)  =  4.4  m,  ay  =  4 1 . 1  pm,  az  =  212.7  pm,  ay'  =  9.7  prad,  and  o2' 
=  18.9  prad..  We  considered  here  the  effect  of  aberrations  only  and  put  aside  any  contribution  from  higher  order  overlapping 
and  surface  scattering.  The  results  will  be  given  in  the  following  subsections  B  and  C. 


B.  Resolving  power 
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To  gain  a  better  understanding  of  the 
resolving  power,  we  have  examined  the 
behavior  of  .iff  1 )  and  M2)  as  a  function 
of  the  ruled  width  W.  This  is  because  the 
resolving  power  is  influenced  directly  by 
the  amount  of  aberrations  in  spectral 
images,  which  increase  with  increasing  W. 

We  constructed  spot  diagrams  and  line 
profiles  for  the  three  gratings  HG(S+S), 

HG(A+S),  and  RG(VLS)  as  a  function  of 
W.  The  ruled  width  was  varied  by 
masking  down  the  ruled  area  stepwise 
from  fV=60  mm  (full  ruled  width)  to 
W=1 .5  mm  by  a  factor  of  0.5,  while 
keeping  the  groove  length  unmasked 
(Z>30  mm).  To  take  into  account  these 
four  different  masked  apertures  on  the 
grating,  we  also  computed  the 
transmittance  of  rays  through  the 
respective  apertures.  For  this,  rays  were 
generated  randomly  in  the  central  cone  of 
the  undulator  until  the  number  of  rays 
diffracted  from  the  unmasked  ruled-area 
reached  1000.  Then,  the  transmittance 
was  given  by  the  ratio  of  1 000  to  the  total 
number  of  rays  generated. 

Figure  8  shows  spot  diagrams  and 
line  profiles  thus  obtained  for  rays  of  X =5 
nm  in  the  cases  of  (a)  HG(S+S),  (b) 

HG(A+S),  and  (c)  RG(VLS).  The  values 
of  transmittance  were  found  to  be  the 
same  in  all  the  cases  (a),  (b),  and  (c); 

28.83%,  55.25%,  88.42%,  and  99.80% 
for  W  =7.5,  15,  30,  and  60  mm, 
respectively.  These  values  were  taken 
into  consideration  in  drawing  the  line 
profiles.  The  full  width,  A  WF,  of  the 
individual  line  profiles  at  their  (4/u2)- 
maximum  points  is  indicated  by  arrows 
and  the  effective  Gaussian  line  profiles  are 
shown  by  dotted  lines.  Before  examining 
the  behavior  of  #(1)  and  M2),  we  should 
estimate  their  statistical  dispersion  arising 
from  random  generation  of  rays.  For  this  we  constructed  spot  diagrams  and  line  profiles  by  tracing  1000  diffracted  rays  for  ten 
different  sets  of  randomly  generated  rays  of  5  nm  through  the  individual  gratings  of  Fig.  8  as  a  function  of  the  width  W  of  the 
unmasked  portion  of  the  ruled  area.  The  resulted  spot  diagrams  and  line  profiles  were  found  to  be  veiy  similar  to  the 
corresponding  ones  in  Fig.  8.  The  mean  values  of  <5*(1)  and  M2)  obtained  are  summarized  in  Table  1.  The  statistical  dispersion 
of  M 1 )  is  much  larger  than  that  of  MX),  but  not  large  enough  to  alter  any  conclusion  to  be  drawn  from  their  mean  values. 
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Fig.  8.  Spot  diagrams  and  line  profiles  constructed  as  a  function  of  the  ruled 
width  W for  the  M-G  type  VLS-PGM  equipped  with  (a)  HG(S+S),  (b)  HG(A+S) 
or  (c)  RG(VLS). 
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Table  I.  Mean  values  of  the  resolving  powers  91(1)  and  91(2)  of  the 
M-G  type  VLS-PGM  as  a  function  of  the  width  W  of  the  unmasked 
portion  of  the  ruled  area. 


For  HG(S+S)  we  notice  in  Table  I  that  with  an 
increase  in  W,  the  gradual  decrease  in  «5¥(1)  forms  a 
sharp  contrast  to  the  rapid  decrease  in  <5¥(2).  We  also 
note  in  Fig.  8a  (1)  that  the  line  profile  rapidly 
develops  shading  (coma-type  aberrations)  with 
increasing  W,  while  A  wF  remains  nearly  unchanged 
and  (2)  that  the  effective  Gaussian  line  profile  nearly 
coincides  with  the  ray-traced  line  profile  for  W  =  7.5 
and  15  mm,  whereas 

those  for  W  =  30  and  60  mm  are  much  wider  than  the 
corresponding  ray-traced  line  profiles. 

HG(A+S)  and  RG(VLS)  exhibit  a  striking 
contrast  to  HG(S+S)  in  their  behavior  with  changes  in 
W.  Regardless  of  the  width  W  of  the  unmasked 
portion  of  the  ruled  area,  both  HG(A+S)  and 

RG(VLS)  provide  nearly  the  same  value  of  -29000  for  M  1 )  and  M2)  and  also  maintain  sharp  Gaussian  type  line  profiles. 
These  findings  suggest  the  need  for  careful  examination  of  the  resolving  power  in  relation  with  the  spectral  purity  of  the  beam 
emerging  from  the  exit  slit. 
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C.  Spectral  purity 


The  spectral  purity  of  the  beam 
emerging  from  the  exit  slit  of  the  M-G 
type  VLS-PGM  was  examined  by 
analyzing  the  frequencies  at  which  rays  of 
various  wavelength  components  passed 
through  the  exit  slit  of  5-pm  wide.  For 
this,  we  assumed  that  the  monochromator 
was  timed  to  a  wavelength  of  5  nm  and 
that  a  continuum  band  which  extends 
over  a  range  of  4.995  nm  <  X  <  5.005  nm 
was  incident  on  the  grating.  Rays  of 
various  wavelengths  in  the  continuum 
band  were  generated  randomly  in  the 
central  cone  of  the  undulator  and  were 
traced  through  the  monochromator 
equipped  with  HG(S+S),  HG(A+S)  or 
RG(VLS)  until  the  number  of  rays  passed 
through  the  exit  slit  reached  1000.  Each 
through  ray  is  identified  by  its  wavelength 
X  and  the  Y  coordinate  (see  Fig.  7)  of  its 
intersection  with  the  plane  of  the  exit  slit. 
Plotting  ( X,Y)'s  for  the  through  rays,  we 
obtain  a  diagram  similar  to  the  spot 
diagram,  which  we  call  the  spectral  purity 
plot  diagram.  The  frequency  distribution 
v(A.)  of  wavelength  components  of  the 
continuum  band  is  estimated  for  the 
through  rays  by  counting  the  number  of 
points  ( X,Y)'s  falling  onto  individual 
parallel  zones  of  5*  10'5  nm  wide  and  5- 
pm  high  in  the  spectral  purity  plot 
diagram.  The  procedure  of  constructing 
the  spectral  purity  profile  is  similar  to  that 
of  making  the  line  profile  from  the  spot 
diagram. 

Figures  9a,  9b,  and  9c  show  the 
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Fig.  9.  Spectral  purity  plot  diagram  and  spectral  purity  profiles  constructed  for 
the  M-G  type  VLS-PGM  as  a  function  of  the  width  W  of  the  unmasked  portion 
of  the  ruled  area. 
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spectral  purity  plot  diagrams 
and  the  corresponding  spectral 
purity  profiles  obtained  for  W 
=  7.5,  15,  30,  and  60  mm  in 
the  cases  of  HG(S+S), 

HG(A+S),  and  RG(VLS), 
respectively.  The  two 
horizontal  dotted  lines  in  each 
spectral  purity  plot  diagram 
indicate  a  slit  opening  of  5 
pm.  The  percentage  T  of  all 
the  rays  generated  that  were 
through-rays  is  indicated  at  the 
top  of  the  respective  spectral 
purity  plot  diagrams.  The 
standard  deviation  ov  of  the 

distribution  v(A)  is  given  at  the  bottom  of  each  diagram.  The  values  of  o ,  are  a  measure  of  the  spectral  purity.  In  Fig.  9,  it  is 
interesting  to  note  the  following  facts.  When  o  v  is  greater  than  1 .7  x  1  O'4  nm,  the  wavelength  band-pass  width  corresponding 
to  a  slit  opening  of  5-pm,  the  spectral  purity  is  degraded  and  the  values  of  M  1 )  and  M2)  are  greatly  different.  In  case  (a)  with 
W- 30  and  60  mm,  ov'sare  larger  than  1.7x10  4  nm,  and  $(1)  =  16100  (or  14300)  and  M2)  ~  6400  (or  2900)  for  W=  30  mm 
(or  60  mm),  respectively.  In  all  other  cases,  ov<1.7xl0~4  nm,  and  $(1)  ~  M2)  regardless  of  the  values  of  W. 

The  spectral  purity  of  the  beam  emerging  from  the  exit  slit  may  also  be  estimated  from  the  ratio  of  the  area  A  m  under  the  main 
peak  between  (5  -  A  A/2)  nm  and  (5  +  A  A/2)  nm  and  the  total  area  Ar  under  the  whole  frequency  distribution  curve,  where  A  A 
is  either  the  full  width  of  the  ray-traced  line  profile  at  its  (4/it  2)-maximum  point  or  2.6423  o  x  according  to  whether  the  resolving 
power  is  defined  by  M\)  or  M2).  Accordingly,  we  differentiate  the  area  A  x  specified  by  M 1 )  from  the  one  designated  by  M2) 
and  denote  the  former  by  1)  and  the  latter  by  AJ. 2).  The  ratios  A^\)IAT  and  A^(2)/AT  were  computed  for  HG(S+S), 
HG(A+S),  and  RG(VLS)  as  a  function  of  IF,  the  width  of  the  unmasked  portion  of  the  ruled  area.  The  results  are  listed  in  Table 
II.  The  table  clearly  shows  ( 1 )  that  in  the  case  of  HG(S+S)  the  ratio  A^\)!AT  rapidly  decreases  with  an  increase  in  IF,  whereas 
AJ2)/At remains  at  ~80%  regardless  of  W  and  (2)  that  in  the  cases  of  HG(A+S)  and  RG(VLS)  all  the  ratios  listed  in  the  table 
remain  nearly  unchanged  as  IF  is  varied  andA^\)IATsn<\Asl2)IAT  are  practically  the  same  for  all  values  of  W. 

The  results  shown  in  Figs.  8  and  9  and  Tables  I  and  II  indicate  that  the  definition  of  M2)  is  more  realistic  as  compared  with 
others  especially  when  the  monochromator  under  consideration  yields  spectral  images  with  asymmetric  profiles. 


Table  II.  Degree  of  the  spectral  purity  of  the  beam  emerging  from  the  M-G  type  VLS-PGM 
as  a  function  of  the  width  W  of  the  unmasked  portion  of  the  ruled  area.. 


IF 

(mm) 

HG(S+S) 

HG(A+S) 

RG(VLS) 

AtfXMr 

A£2)!At 

A^\)IAr 

AJ2)/AT 

A^I)/AT 

A^2)/At 

7.5 

84% 

84% 

85% 

85% 

86% 

85% 

15 

73 

76 

86 

85 

85 

84 

30 

58 

74 

85 

84 

86 

85 

60 

49 

79 

84 

85 

82 

85 

7. 3.  Thermal  deformation 

The  heat  load  exerted  on  beamline  optics  by  undulator  radiation  can  distort  their  surfaces.  The  effect  is  the  greatest  on  the 
first  optical  element,  which  is  usually  a  mirror.  The  performance  of  such  a  thermally  deformed  mirror  can  be  simulated  by  means 
of  ray  tracing.  The  thermal  deformation  is  calculated  in  four  steps:  (1 )  calculation  of  individual  harmonic  power  densities  of 
undulator  radiation  incident  on  the  mirror  surface  by  means  of  a  code  such  as  URGENT  [8  ],  (2)  calculation  of  the  harmonic 
power  densities  absorbed  at  each  mesh  element  on  the  mirror  surface,  (3)  calculation  of  thermal  deformation  by  means  of  a 
finite-element  code  (e.g.,  ANSYS,  a  commercial  soft  ware),  and  (4)  derivation  of  an  analytic  expression  of  A  £TO  by  fitting  the 
surface  deformation  data  with  a  polynomial  of  the  form 

-  £  Aijw  '/y>  0  *  i*J  *  6>  (20) 


where  A  2s  are  the  deformation  coefficients. 

As  an  example  we  show  the  result  of  a  simulation  made  on  the  thermal  deformation  of  a  spherical  mirror  exposed  to 
undulator  radiation.  The  machine  parameters  assumed  for  the  undulator  are:  wavelength  of  the  fundamental  radiation  =  4.5  nm; 
undulator  period  length  Au  =  39  mm;  number  of  magnet  period  N  =  123;  deflection  parameter  K =1.4057;  electron  energy  /  me 
c 2  (me,  electron  mass;  c,  velocity  of  light)  =  y  =  2395.08;  rms  transverse  size  of  the  electron  beam  in  tire  horizontal  direction 
oy=  0.330  mm;  rms  transverse  size  of  the  electron  beam  in  the  vertical  direction  oz  =  0.063  mm;  rms  angular  divergence  of  the 
electron  beam  in  the  horizontal  direction  ay'  =  30  prad;  rms  angular  divergence  of  the  electron  beam  in  the  vertical  direction 
o/  =  16  prad.  The  mirror  is  assumed  to  accept  undulator  radiation  of  0.5  mrad  (h)  x  0.5  mrad  (v)  through  a  beam  defining 
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aperture.  It  is  also  assumed  that  the  mirror  is  a  directly  water-cooled  metallic  spherical  concave  mirror  having  the  following 
specifications:  radius  of  curvature  R  =  75  m;  dimensions  =  178  mm  (W)  x  66  mm  (H)  x  28  mm  (D);  blank  material  =  brazed 
assembly  of  Glidcop™  and  OFHC  copper,  whose  surface  is  plated  with  electroless  nickel;  and  water-cooling  channel  =  machined 
series  circuit  with  a  rectangular  cross  section  (hydraulic  diameter  =  0.4763  cm)  for  turbulent  water  flow  (flow  rate  =  5.53 
liters/min  and  Reynolds  number  >  4000). 

The  surface  deformation  of  the  mirror  M  is  calculated  by  using  URGENT  and  ANSYS.  URGENT  calculation  shows  that 
the  mirror  receives  a  total  power  of 209  W  and  absorbs  127  W  on  its  surface.  The  results  of  ANSYS  calculation  are  curve  fitted 
to  Eq.  (20).  This  yields  the  thermally  deformed  surface  of  the  spherical  mirror  M  as 

i  .  R  -  R[1  -  (w2*  l2)/R2]m*  3.09499  xlO-4-  3.39872  xl(T8w2-  8.21102  xHT8/2 

*  4.75112  x1(TI2h’4+  4.833510  xlO'11*2/2-  4.880480  xl0“7/4-  3.812470  xl0"16w  6  'irP}  (21) 

-  1.569350  xlO'14*4/2*  5.484390  xio ~nw2lA.  2.030340  xl(T8/6, 

where  R  is  the  radius  of  curvature  of  the  mirror.  '  J  ;V>'!  ‘ 

To  evaluate  the  performance  of  this  deformed  mirror  by  means  of  ray-traced  spot  diagrams,  it  is  required  to  incorporate  Eqs. 
(15)  and  (2 1 )  into  an  exact  ray  tracing  procedure.  The  ray-tracing  procedure  can  easily  be  modified  to  meet  the  requirement 
by  following  the  method  given  in  Ref.  9.  Rays  from  the  central  cone  of  the  undulator  can  be  generated  by  a  simulation  code  [10] 
that  takes  into  account  the  machine  parameters  and  emission  probability.  '  ’l  1  u'*'i 

8.  Demagnifying  Projection  Optics  •  m  tajnmw j 

Examples  are  given  for  the  tolerance  for  the  figure  errors  in  the  individual  components  of  projection  optics,  the  effect  of 
figure  errors  in  projection  optics  on  the  resolution  of  mask  pattern  images,  and  simulation  of  Foucault  knife  edge  testing  for 
assembly  errors  [1 1].  "  -f*f|  A 

8.  1.  Tolerance  for  the  figure  errors  of  ellipsoidal  mirrors 

A  measure  of  figure  errors  that  are  allowable  to  achieve  near  diffraction  limited  performance  of  a  mirror  system  is  provided 
by  Rayleigh's  quarter  wavelength  rule  [12]  and  the  Marechal  condition  [13].  These  tolerance  criteria  give  A/8  peak-to-valley 
(P-V)  and  A/28  rms  for  each  mirror  of  the  two-mirror  system:  1 .6  nm  P- V  and  0.46  nm  rms  for  A  =  1 3  nm.  These  tolerances 
often  are  too  tight  in  practice.  It  is  therefore  desirable  to  develop  a  means  of  obtaining  a  practically  reasonable  tolerance  for  the 
figure  error.  One  approach  to  this  task  would  be  the  use  of  ray  tracing.  We  have  tried  to  utilize  the  method  described  in  Section 
6  for  exploring  the  possibility.  For  this  we  used  the  NTT's  projection  optics  as  an  example. 

The  optical  configuration  of 
the  NTT's  EUV  projection  litho¬ 
graphy  system  is  illustrated  in  Fig. 

10.  The  system  consists  of  the 
illumination  optics  for  producing 
a  ring  field  on  a  reflection  mask 
and  the  projection  optics  for 
imaging  a  reduced  pattern  of  the 
mask  on  a  wafer. 


Synchrotron  radiation  of  40 
mrad  (h)  x  4  mrad  (v)  from  a  600 
Mev,  500  mA  (1 .2  A  max)  super¬ 
conducting  compact  electron 
storage  ring  at  the  NTT  SR  facility 
is  shaped  to  a  parallel  beam  of  20 
mm  (h)  x  10  mm  (v)  by  a  concave 
toroidal  mirror  and  a  convex 
toroidal  mirror  and  then  focused 
on  the  reflection  mask,  RMK,  by  a 
rotatable  concave  cylindrical 
mirror  and  a  concave  toroidal  mirror. 
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Fig.  10.  Optical  configuration  of  the  NTT's  EUV  projection  lithography  system. 

The  focused  beam  illuminates  a  ring  field  of  5  mm  (h)  x  1  mm  (v)  on  RMK,  the  radius 
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of  the  ring  image  being  62.5  mm.  The  ring  field  is  expanded  to  100  mm  (h)  x  1  mm  (v)  by  oscillating  the  cylindrical  mirror 
through  angles  of  ±53. 13°,  and  it  is  further  extended  to  100  mm  (h)  x  1 25  mm  (v)  by  translating  RMK  vertically. 

The  illuminated  area  of  RMK  is  projected  on  the  wafer  by  a  demagnifying  optics  with  a  numerical  aperture  of  0. 1  and  a 
magnification  of  1/5.  The  demagnifying  optics  is  of  telecentric  and  consists  of  a  convex  ellipsoidal  mirror  Ml  of  k  =  ~10  and 
a  concave  ellipsoidal  mirror  M2  of  k  =  ~0. 1 ,  k  being  the  conic  coefficient  (see  Fig.  10).  This  coaxially  symmetric  two-mirror 
system  was  chosen  because  of  its  simplicity,  reasonable  aberration  balancing,  high  throughput,  relatively  large  fabrication 
tolerance,  and  easy  assembling.  The  diameters  of  Ml  and  M2  are  60  mm  and  1 50  mm,  respectively,  and  both  the  mirrors  have 
almost  the  same  radius  of  curvature,  —500  mm.  The  surfaces  of  Ml  and  M2  are  very  close  to  sphere:  the  maximum  deviation 
of  the  surface  figure  of  Ml  (or  M2)  from  its  reference  spherical  surface  is  less  than  0.6  (or  2.0)  pm.  An  addition  of  vertical 
translation  of  the  wafer  to  the  motions  of  the  rotatable  cylindrical  mirror  and  RMK  will  produce  the  image  of  20  mm  (h)  x  25 
mm  (v)  on  the  wafer. 


°se(M2)  -  0 

°se(M  1 )  =  0.05  prad 


aSE(M2)  =  0 
aSE(Ml)  =  0.1  prad 


cjse(M2)  =  0 
ose(M1)  =  0.2  prad 


c 

=L 


pm  /  division 


Fig.  1 1.  Ray-traced  images  of  the  source  pattern  C  on  the  mask.  The  source  pattern  is  of 
Fig.  4  with  a  =  0.5  pm.  oSE  is  the  rms  slope  error  of  the  mirror  specified  in  parentheses. 
Ml  and  M2  are  the  convex  and  the  concave  ellipsoidal  mirror. 


To  estimate  a  practically 
reasonable  tolerance  for  the 
figure  error  of  the  demagnifying 
projection  optics,  we  follow  the 
procedures  described  in 
Sections  4.2  and  6.  We  consider 
a  0.5-pm  line-and-space  grid 
pattern  on  the  reflection  mask 
and  assume  that  the  space 
portions  of  this  grid  pattern  emit 
light  into  the  entrance  pupil  of 
the  system.  We  choose  three  3 
x  3  pm2  areas.  A,  B,  and  C.  The 
areas  A  and  B  are  taken  at 
around  the  two  comer  points  at 
the  right  (or  left)  of  the  ring  field 
of  5  mm  (h)  x  1  mm  (v)  formed 
on  RMK  and  the  area  C  is  taken 

at  around  the  center  of  the  ring  field.  We  consider  these  area  as  distributed  light  sources ,  each  containing  9  elemental  sources 
of  0.5  pm  x  0.5  pm.  Ten  thousand  rays  generated  from  each  one  of  these  three  sources  A,  B,  and  C  were  traced  through  the 
convex  ellipsoidal  mirror  Ml  and  the  concave  ellipsoidal  mirror  M2.  The  spot  diagrams  thus  obtained  with  the  three  sources 
were  found  to  be  indistinguishable  from  one  another  as  far 
as  the  same  slope  errors  are  assumed  for  the  projection 
system.  In  Fig.  11  are  shown,  as  an  example,  the  spot 
diagrams  constructed  for  the  source  C  and  for  certain 
values  of  slope  errors  assumed  on  Ml  and  M2. 

While  constructing  such  spot  diagrams  with  assumed 
slope  errors  of  Ml  and  M2,  we  can  calculate  from  Eq.  (15) 
the  corresponding  figure  errors  of  Ml  and  M2.  Figure  12 
shows  the  relation  between  the  rms  slope  error  and  the  rms 
figure  error  thus  obtained  for  Ml  and  M2  of  the  NTT's 
projection  optics.  It  is  seen  that  the  rms  figure  error  is 
affected  by  the  rms  slope  error  more  strongly  for  M2  than 
for  Ml. 

To  estimate  the  tolerance  for  the  figure  error,  we  need 
information  on  the  intensity  distributions  of  the  images  in 
the  horizontal  and  vertical  directions.  Such  information  can 
be  obtained  from  the  spot  diagrams  by  constructing  image 
pattern  profiles  in  the  two  directions.  The  image  pattern 
profiles  were  constructed  by  counting  the  number  of  spots  j-^g  glope  error  (|irad) 

falling  into  individual  vertical  or  horizontal  zones  of  0.006- 

pm  in  width  in  the  image  plane.  The  intensity  profiles  of  Fig  j  2  Relation  between  the  mis  slope  error  and  the  mls 

spot  diagrams  constructed  for  Ml  and  M2  with  nine  sets  of  flgure  error  m  Ule  case  of  ^  N1Ts  project,on  system 

figure  errors  are  shown  in  F,g.  13.  consisting  of  the  ellipsoidal  mnrors  Ml  andM2. 
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Fig.  13.  Spot  diagrams  and  their  intensity  profiles  constructed  for  the  demagnifying  projection  system  with  various  sets  of  the 
rms  figure  errors  of  Ml  and  M2.  '  |||| 

In  Fig.  1 3  the  rms  figure  errors  are  assumed  to  be  ore  =  0,  1 .0,  and  1 .8  nm  for  the  convex  ellipsoidal  mirror  Ml  and  o 
=  0, 1.0,  and  1.5  for  the  concave  ellipsoidal  mirror  M2.  The  ">45%"  criterion  was  applied  to  the  modulations  estimated  from 
Fig.  13  and  other  similar  spot  diagrams.  The  result  shows  that  any  combination  of  ofe(M1)>  1.2  nm  and  oFE(M2)  >  1.5  nm 
yields  a  modulation  less  than  45  %.  We  therefore  adopt  1 .2  nm  and  1 .5  nm  as  the  tolerances  for  the  rms  figure  errors  of  Ml  and 
M2,  respectively. 


The  surface  figures  of  Ml  and  M2  fabricated  were  measured  with  a  Zygo  Mark  IV  interferometer.  The  rms  (or  Peak-to- 
Valley,  P- V)  values  of  the  figure  errors  were  found  to  be  ~1 .8  nm  (or  ~5.5  nm)  and  ~1 .5  nm  (or  ~4.4  nm)  for  Ml  and  M2, 
respectively.  The  spot  diagram  and  its  intensity  profiles  constructed  with  these  tolerance  values  are  shown  at  the  lower  right 
of  Fig.  1 3 .  The  figure  clearly  shows  that  a  0. 1  -pm  line-and-space  pattern  cannot  be  resolved  with  the  projection  optics 
fabricated,  in  agreement  with  the  experimental  result.  This  suggests  the  need  for  improved  surface  figures  of  Ml  and  M2.  As 
stated  above,  the  present  simulation  gave  the  tolerances  of  1 .2  nm  (or  0. 1 1  prad)  and  1 .5  nm  (or  0.04  prad)  for  the  rms  figure 
error  (or  rms  slope  error)  of  Ml  and  M2,  respectively.  These  values  would  ease  the  Mareshal  condition  of  0.46  nm  rms  for 
both  Ml  and  M2,  but  their  reliability  still  remains  to  be  proved  experimentally. 


8.  2.  Simulation  of  Foucault  knife  edge  testing  for  assembly  errors 

This  subject  is  reported  by  T.  Haga  and  M.  C.  K.  Tinone  in  the  Proceedings  of  this  workshop.  The  interested  reader  is 
referred  to  their  article  entitled  "Numerical  analysis  of  EUV  optics  assembly  using  at-wavelength  Foucault  testing"  in  this 
volume. 
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